Superdiffusivity of quantum walks: A Feynman sum-over-paths description 
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Quantum walks constitute important tools in different applications, especially in quantum algo- 
rithms. To a great extent their usefulness is due to unusual diffusive features, allowing much faster 
spreading than their classical counterparts. Such behavior, although frequently credited to intrinsic 
quantum interference, usually is not completely characterized. Using a recently developed Green's 
function approach [Phys. Rev. A 84, 042343 (2011)], here it is described - in a rather general 
way - the problem dynamics in terms of a true sum over paths history a la Feynman. It allows 
one to explicit identify interference effects and also to explain the emergence of superdiffusivity. 
The present analysis has the potential to help in designing quantum walks with distinct transport 
properties. 

PACS numbers: 03.67.Lx, 05.40.Fb 



I. INTRODUCTION 

Quantum walks (QWs), a quantum version of classical 
random walks (CRWs) is a relatively simple class 
of systems, yet containing almost all the essential as- 
pects of quantum mechanics @, H| . They can be used 
to model a large number of phenomena [J, such as: 
energy transport in biological systems [(|; Bose-Einstein 
condensates redistribution Q ; quantum phase transition 
in optical lattices Q; and decoherence processes 0. But, 
certainly where QWs have attracted more interest is in 
quantum computing [l(| EH- In fact, QWs allow the 
development of new quantum algorithms (12j , which of- 
ten display much better performance than their classical 
siblings [12, [H| . 

QWs usefulness in applications is in great part due to 
their unusual transport properties. For instance, they 
present exponentially faster hitting times [F3 - [la ] (the 
time necessary to visit any vertex in the system graph 
space), an im por tant feature for searching in discrete 
databases [l7l Ilcf. Such faster spreading compared to 
CRWs is usually attributed to interference [1, HU , a 
key ingredient in implementations (20 [ [2 lT| and believed 
central to explain distinct behaviors [22|. However, ex- 
actly how these effects emerge in QWs usually is not to- 
tally characterized [23j . posing challenges as to how one 
could properly link the high degree of entanglement in 
QWs [24j ] with interference. Furthermore, since interfer- 
ence actually comes from a high proliferation of paths 
(after all, QWs are associated to the idea of CRWs 25] ) , 
it also bears on the problem of how decoherence [26| can 
make the "quantum trajectories" to become classical [27| . 

In trying to understand interference in QWs, a path 
integral-like treatment would be appropriate. Actually, 
a few interesting works along this line have been proposed 



28]. However, they address the problem from a different 
perspective^using combinatorial analysis to compute fi- 
nal states [291 ] , but not considering intermediary steps in 
terms of Feynman's history of trajectories [301 ] . Hence, 
interference is not made truly explicit. 

In the present work we show how quantum interference 
determines QWs uncommon diffusive properties. To this 
end, the exact Green's function [3l[ - given as a general 
sum of paths - is written in a closed analytical form. 
Then, we describe how to calculate relevant quantities in 
a way identifying the trajectories superposition contri- 
butions. To concretely illustrate the approach, we show 
that the usually observed: (a) complicated oscillatory 
behavior of the probability distribution for visits at dif- 
ferent sites; and (b) the process dispersion dependence on 
time; are associated to the complex multiple reflections 
and transmission patterns of the system evolved paths. 

Finally, we mention some important technical aspects. 
There are several ways to formulate QWs, all defined in 
discrete spaces (graphs) [2J. Also, time may be either a 
continuous [14] or a discrete variable. In the latter, the 
major formulations are coined [HJ and scattering [33[ 
QWs. Continuous time and coined are related through 
appropriate limits (34j , whereas coined and scattering are 
unitarily equivalent in any topology and for arbitrary 
transition amplitudes (35| . Hence, we consider only scat- 
tering quantum walks (SQWs), keeping in mind that our 
finds can be extended to such other constructions as well. 
Moreover, the Green's function method considered here 
[3lj is valid for any graph topology. Although for our 
purposes we address QWs on the line, avoiding extra 
and unnecessary mathematical complications, we men- 
tion that the same type of analysis would likewise work 
in more complex networks. 
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II. THE SUM OVER PATHS DESCRIPTION 

We assume an undirected ID lattice of equally spaced 
vertices labeled in Z, Fig. [TJ Pairs of neighbor vertices 
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FIG. 1. QWs graph structure in ID. At each edge there are 
two basis states, e.g., \ + ,j) and |— , j — 1) (schematically rep- 
resented by arrows) for the edge between j — 1 and j. In detail 
the vertex dependent scattering quantum amplitudes. 
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FIG. 2. For 5, the three possible situations for the relative 
positions of the initial and final, i and /, edges. 



are joined by a single edge. To each edge we ascribe two 
basis states. For instance, in Fig. [T]for the edge between 

I 



j — 1 and j (j and j -I- 1) we have | + 1, j) and | — l,j — 1) 
(| + 1, j + 1) and | - 1, j)). Therefore, the full set {\a,j)\ 
spans all the possible system states \tp). The quantum 
numbers a = ±1 (hereafter for short ±) represent the 
propagation direction along the lattice (or graph). The 
discrete dynamics is given by the one step time evolution 
operator U, such that the state at times m + 1 and m 
are related by \ip(m + 1)) = U\tp(m)}. For an arbitrary 
phase z = exp[^7], we have [3l[ 



U\a, j) = z (tf \a, j + a) + r> } \-a,j-a}). (1) 
Here [H (with < r^tj < I and < ^ < 27r ) 



tS ±} = 
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rj exp[«0;./J 



±7T, 



(2) 



guarantee the evolution unitarity. The r's and i's can 
be understood as the vertices reflection and transmission 
quantum amplitudes (see Fig. [1]). 

The problem is fully described by the Green's function 
approach in [3l| . Consider the walk starting at the edge 
i (between the vertices j — 1 and j) with the state di- 
rection a and finally getting to the edge / (between the 
vertices j ± n and j ± (n — 1), for any n — 0, 1, 2, . . .). 
The three possible situations are illustrated in Fig. 2. 
Then, the most general exact expression for the Green's 
function, representing the transmission probability am- 
plitude {a, i} —> f reads (s = — 1 (+1) [in short =f] for / 
to the right (left) of i and s = for / equal to i) 
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The composed coefficients and T-^, functions of 

the individual amplitudes r~- s and ^ s > are obtained 
from the following recurrence relations [36j (/i- = j — 
(s + l)(n- l)/2 and /z+ =j - 1 - (s - l)n/2 for s ^ 0) 
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In Eq. <j3j) it is not specified what is the final direction 
quantum number, v, when arriving at /. In fact, it in- 
cludes both cases once v = a (—a) corresponds to the 
term I (zR) in the second (. . .) in the numerator of Eq. 



([3]). There are different contexts for which we may seek 
the amplitude transition {er, i} — > f. Common ones are: 
(i) exactly after m — M time steps; and (ii) when the 
system never visits vertices further to the left and to the 
right than, respectively, j = Ji and j = J r (e.g., for first 
passage time calculations) . In both we just need two ex- 



?(-) _ „(-) 



tra relations for Eq. Ry> = ry and Rj 
Moreover, for (i) we have J; = (j— 1) — [(M+sn— Si s )/2] 
and J r = j+[(M—sn+Si s )/2], with [x] the integer part of 
x and n taken consistently. Finally, for Eq. Q obviously 
|V(0)) = |a,i + (<r-l)/2). For |^(0)) =E^jW,j), the 
correct Green's function would be G = J2 c °,j Gf,{<y;i\& •}■ 
The above exact expression is derived from a sum 
over infinite many "scattering paths" jU starting 
and ending at the edges i and /. Its advantage is that 
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all the possible quantum walk trajectories are "com- 
pacted" into a closed formula. So, distinct interference 
phenomena can be extract from Q. Indeed, as demon- 
strated in 31], this is achieved in a rather systematic 
way by means of two differential operators. The proba- 
bility for {a, i} — > f in exactly m time steps is given by 

P{a,i}^f{m) = \S m Gf,{a,i}\ 2 , With S„ 



m=0 m=l 



1 d m 

ml dz n 



, = n the 



Step Operator. To see it, we note that [3l| S m Gf { CT %} 



.s . p 



V s . p ., for each V s . p . being the contribution of a tra- 
jectory from i to / in m steps (38|. Interference comes 
into play when we take the modulus square of such ex- 
pansion. Also, any specific V follows from PpG, for the 
Path Operator (superscripts ± omitted for clarity) 



Pv = 



TT 

- LJ - mt! mil 

k.iev 



gm k gmi 

dt k k dr T l 



(5) 



In Eq. ([5]) the t^s and r;'s are the scattering amplitudes 
(appearing m^ and mi times) characterizing the path as- 
sociated to V [311 ] . 

With this mathematical 'machinery', below we can 
make an analysis of sum over paths for QWs. 



III. RESULTS AND DISCUSSION 



Suppose the initial state \ijj(0)) = so we write 




j-5 j-4 j-3 j-2 j-1 j j+1 j+2 j+3 j+4 j+5 



FIG. 3. For each m, an arrow — > (<—) pointing to the vertex 
j' indicates that the basis state | + , j'} (|— is present in 
the expression for \ip(m)) — U m |+, j). The number of arrows 
of a given type equals the number of distinct paths leading to 
the corresponding \u,j'}. 



(first step), and finally goes to j — 2 (second step), now 
heading left. Note that in terms of a quantum scattering 
process, it represents a reflection from the vertex j (rj- ) 
and then a transmission through the vertex j — 1 (t^Ci). 

By applying S m on G{v,f},{a,j} an( i afterwards simply 
setting = tj =1 Vj, one directly finds that the 
number of paths leading to {cr,j} — > \v, j'} after m steps 
is given by the binomial coefficient [39( 



U m \m) = Mm)) = 



j =j+m 

" E E 

j'=j-m v=± 



(6) 



Above, some a's are zero since certain basis states are 
absent, e.g., \—,j + m—l) cannot be reached in m steps. 
In fact, exactly 2m a's are not null. 

Next, to make contact with CRWs, we observe that by 
leaving from the edge corresponding to {a, j} there is a 
determined number of trajectories (eventually none) fi- 
nally getting to specific edges in exact m steps. Thus, 
the total number of paths ending up in any possible 
\v,j') is 2 m . Since p — |a| 2 , the a's can be given as the 
sum of the quantum amplitudes 38] of all paths yielding 
{a,j}^{v,f}. 

As an simple example, for |^>(0)) = \+,j) the Fig. 
[3] schematically illustrates the basis states expansion of 
\tp(m)) (up to m = 5). Consider m = 2, then 

z - 2 \m) = z- 2 u 2 \+,j) 

= a-,j-2 - 2) +a +ij 

+a- <j \-,j)+a +<j+2 \+,j + 2). (7) 

The presence of, say, \—,j — 2) in the expression for 
|V'(2)), Eq. ([7]), is represented in the m = 2 case of 
Fig. Oby an arrow pointing to the vertex j — 2. More- 
over, its number, here just a single arrow, means there is 
only one path getting to — 2) from if m = 2: 

a trajectory which initially heading right at j (since 
| V(0)) = \+,j)), reverses its direction at j, goes to j — 1 



m — 1 

m+j'-j _ r 



(8) 



As a simple check, one can test Eq. (j8|) with the schemat- 
ics in Fig[3l Moreover, from the mapping between SQWs 
and coined QWs in [35|, the number of paths to a given 
j' state for the latter QW formulation is trivially derived 
from N v .y as 



'coined 



E^: 



m+j'-j 
2 



m-j'+j 
2 



(9) 



which agrees with the formula in Ref. [4(| (with j = 0). 

Assume any path taking, regardless the order, 
(dM) steps to the left (right). It would lead the system 
to \v,j' = j + d( + ) — S^). Reversing this reasoning, 
consider a fixed m = d' + ' > \Aj\, with Aj = j'—j. 

Paths for which eft ' = (m± Aj)/2 are both integers will 
result in j —> j' . To obtain all such paths, we should 
consider Q for J; = j — and J r = j + S + \ The 
contribution from each path to a given coefficient in Eq. 
([6]) will involve exactly m position dependent amplitudes 
Tj's and t^'s. In this way, the actual procedure to 
calculate the a's is to compute a u ,j' = SmG{u,j'},{(r,j\i f° r 
Ji, J r , d^' as above. 

For complete arbitrary rj's and tj s and for m in 
the hundreds, any available computer algebra system can 
be used to obtain the a's as explained. Actually, vertex- 
dependent quantum amplitudes can give rise to a great 
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diversity of diffusive properties [4i(. So, the present pro- 
cedure may be useful to test distinct QWs models trans- 
port features, helping to choose sets of rj- 's and 4 s 
more appropriate in different applications (examples to 
appear elsewhere). 

However, a real surprise is for the situation when su- 
pcrdiffusion takes place even for j independent quantum 
amplitudes and when at each single step the QW resem- 
bles an unbiased classical walk (i.e., 50-50% probability 
to go right-left) 2}. In the following, we show how in- 
terference can fully explain this apparently non-intuitive 
behavior. 

= and from the above pre- 
min{d^ - S av , - I}) 



T7 (±) 

_bor Tj 



scription, we get {n sup 

fn C n , 



n—n su 
n= — 8 (T 



fn 



d(- CT > - 1 



C n 



x r t (-ir)]iJ(-'')-n-lr r W]n+l i 
Furthermore, using Eq. ([2]) 



exp[z0] t m ( 7 ) (-1)" 



(10) 



(11) 



with 



a global phase (unimportant here) which depends 
on j, j', a, is, to and tj>^ . In Eq. (fit))) . /„ gives the 
number of distinct paths yielding a same amplitude C n 
to the a's. This is possible because different paths cor- 
respond to a different order of scattering processes along 
the lattice. Nevertheless, if the final set of scattering's co- 
incides, the resulting amplitudes C n are equal. The total 
number of paths for a u ,j' is N v ji = J2n /« = (d(™-<5 CT J' 
which agrees with Eq. ([8]). 

Particularly important in Eq. (fTTj) is the factor (—1)™, 
arising from the phases difference, Eq. ((2J, between re- 
flections and transmissions in a trajectory. In fact, for 
each path the number of directions change along the way 
is 2n + 1 + 8 av . Therefore, distinct paths may contribute 
with distinct signals (through (—1)") to the sum in Eq. 
(|10j). leading to constructive or destructive interference. 

Lastly, in the "unbiased" case of r = t = l/v2, i.e., 
50%-50% reflection-transmission probability in each ver- 
tex (for a similar coined case see, e.g., 40]), Eq. (TTU)) 
reduces to 

auj, = exp[z0] 2- m / 2 {-2 m 6 mdM + [d^]*" 
x 2 Fi(-dW + S au ,-d^ + l;l + <y w ;-l)}, (12) 

with 2F1 the Gaussian hypergeometric function. To il- 
lustrate this formula, we consider the a's for the final 
states \+,j + 1) and \+,j + 3) in the case to — 5 of Fig. 
[3] (thus b av = I). From Eq. (TT2"]) we get = and 

l a +-j+3| 2 = 1/2- To understand why, note that from /„ 
in Eq. (|I0l) or by inspecting Fig. |31 we find there are 
six (four) possible paths leading to |+, j + 1) + 3)). 

For j' = j + 1, three paths have two (n = 0) direction 
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4. Up to a global phase, the dimensionless coefficient 
a v ji, Eq. ()12l) . as function of the quantum number j' for 
I V(0)) = | + , j), m — 100, (a) v = +, and (b) v = — . Since m 
is even, a Vl j> = if j' — j is odd. The notorious 0| amplitudes 
asymmetry because the particular initial state arises only for 
Situ ~ 1 (case (a)), (c) The linear dependence of the standard 
deviation A on the discrete time m. In the insets \a V ji\ 2 
and the corresponding probabilities for an unbiased classical 
random walk (dashed curves) vs. j' — j. 



changes and three have four (n = 1). The phases are 
then, respectively, (—1)° = 1 and (— l) 1 = —1. Hence, 
these two groups of paths suffer destructive interference. 
On the other hand, for j' = j + 3 there are four pos- 
sible paths, all with two direction changes (n = 0) and 
thus with a same phase. The paths therefore build up a 
relatively high amplitude. 

The above results can also explain two typ ical and im- 
portant behaviors observed in QWs 0, |22| | (see Fig. @] 
(a)-(c)): (i) for usual CRWs, the probabilities for the 
particle location are Gaussian distributed, with a stan- 
dard deviation of y/m (insets of Fig. H](c)). On the other 
hand, quantum mechanically the |a| 2 's, representing the 
particle distribution along the graph, are not spatially 
concentrated; (ii) ay vs. j 1 presents stronger oscillations 
for the j n s far away from the initial j, a pattern usually 
without a classical analog. 

In fact, both (i) and (ii) originate from a similar mech- 
anism. For a fixed large to, the number of trajecto- 
ries N v .y leading to j' is large (small) if |Aj| is small 
(large). In the classical case, since there are no interfer- 
ence, the probabilities are directly proportional to N v j> 
and the Gaussian distribution naturally emerges (recall 
that binomial distributions, c.f. Eq. converge to 

Gaussians). In the quantum case, the many cancella- 
tions coming from opposite signals for distinct groups of 
trajectories, Eq. (TTTj) . prevents the probabilities at |Aj| 
small to be much higher than those at larger |Aj|, Fig. 0] 
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(a)-(b). Hence, a more balanced distribution among the 
states j"s is obtained. By the same token, the smooth 
(strong oscillatory) behavior for |Aj| small (large) is due 
to the fact that varying j' in such interval will propor- 
tionally cause a small (large) change in the number of 
paths contributing to ay. This results in a slow (rapid) 
variation of a' as function of j'. 

Thus, the observed system fast spreading, e.g., in the 
unbiased case characterized by a linear dependence on 
m for the standard deviation (pj> = |a +J '| 2 + |a_.j'| 2 ): 

A = V /E / (J'-J) 2 p j '-(E / (J V -J>i') 2 (Fig- S (c)); 

is due to (i)-(ii). In their turn, (i)-(h) are a direct conse- 
quence of intricate interference effects among paths with 
different phases. 

IV. CONCLUSION 

Summarizing, our contribution here has been twofold. 
First, we propose a distinct approach - based on a true 
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